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Abstract 

A detailed investigation of the scaling properties of the fully finite 0{n) sys- 
tems, under periodic boundary conditions, with long-range interaction, de- 
caying algebraically with the interparticle distance r like r~ d ~ a , below their 
upper critical dimension is presented. The computation of the scaling func- 
tions is done to one loop order in the non-zero modes. The results are ob- 
tained in an expansion of powers of y/e, where e = 2a — d up to 0(e 3 ^ 2 ). 
The thermodynamic functions are found to depend upon the scaling variable 
z = -Rt/~ 1 / 2 L 2_,?_E / 2 , where R and U are the coupling constants of the con- 
structed effective theory, and L is the linear size of the system. Some simple 
universal results are obtained. 
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I. INTRODUCTION 



The theory of continuous phase transitions is based on the hypotheses that at temper- 
atures close to the critical T c , there is only one dominating length scale related with the 
critical behavior of the system. Because of the divergent nature of the correlation length as 
the critical point is approached, the microscopic details of the system becomes irrelevant for 
the critical exponents describing the singular dependence of the thermodynamic functions. 
This intuitive picture is based on the grounds of the renormalization group treatment of 
second order phase transitions. 

Scaling is a central idea in critical phenomena near a continuous phase transition and 
in the field theory when we are interested in the continuum limit [p]]. In both cases we are 
interested in the singular behavior emerging from the overwhelming large number of degrees 
of freedom, corresponding to the original cutoff scale, which need to be integrated out leaving 
behind long-wave length which vary smoothly. Their behavior is controlled by a dynamically 
generated length scale: the correlation length Such a fundamental idea is difficult to test 
theoretically because it requires a study of a huge number of interacting degrees of freedom. 
Experimentally, however, one hopes to be able to study scaling in finite systems near a 
second order phase transition. Namely the system is confined to a finite geometry and the 
finite-size scaling theory is expected to describe the behavior of the system near the bulk 
critical temperature (for a review on the finite-size scaling theory see Ref. [@,|3|]). 

The (9(n)-symmetric vector models are extensively used to explore the finite-size scaling 
theory, using different methods and techniques both analytically and numerically. The 
most thoroughly investigated case is the particular one corresponding to the limit n = oo 
(this limit includes also the mean spherical model) ||. In this limit, these models are 
exactly soluble for arbitrary dimensions and in a general geometry. These investigations 
were devoted exclusively to systems with short (including nearest neighbors) as well as long- 
range forces decaying with the interparticle distance in a power law. For finite n the most 
frequently used analytical method is that of renormalization group However this is 

limited to the case of short-range interaction. The crossover from long to short-range forces 
was discussed in Ref. [|| , where it has been found that renormalized values of the temperature 
and the coupling constant are continuous functions of the parameter controlling the range 
of the interaction, when this approaches the value 2 characterizing the short-range force 
potential. The case of pure long-range interaction was investigated very recently in Ref. 
(a comment on the method and the results obtained there is presented in Section |TV|) . In the 
mean time, a special attention was devoted to the investigation of finite size-scaling for the 
mean-spherical model with long-range interaction (for a review see Ref. |7j and references 
therein) . 

In recent years there has been an increasing interest in the numerical investigation of the 
critical properties of systems with long-range interaction decaying at large distances r by a 
power-law as r~ d ~ u , where d is the space dimensionality and a is the parameter controlling 
the range of the interaction. The mostly used technique for this achievements is the Monte 
Carlo method. This method was used to investigate the critical properties of Heisenberg 
ferromagnetic systems |§ as well as Ising models |9,10[|. Nevertheless all the analysis there 
was concentrated on systems with classical critical behavior in the sense that the critical 
exponents are given by Landau theory. 



2 



In this paper we present a detailed investigation of the finite-size scaling properties of 
the field theoretic 0{n) vector <y? 4 model with long-range interaction. We will also check 
the influence of the interaction range on the critical behavior. These interactions enter the 
exact expressions for the free energy only through their Fourier transform, which leading 
asymptotic is U(q) ~ q a , where a* = min(cr, 2) ||11|| . As it was shown for bulk systems 
by renormalization group arguments a > 2 corresponds to the case of finite (short) range 
interactions, i.e. the universality class then does not depend on a JTTIJT^] . Values satisfying 



< a < 2 correspond to long-range interactions and the critical behavior depends on 
a. With the renormalization group treatment it has been found that the critical behavior 
depend on the small parameter e — 2a — d, where 2cr corresponds to the upper critical 
dimension [|TT | . According to the above reasoning one usually considers the case a > 2 as 



uninteresting for critical effects, even for the finite-size treatments So, here we will 

consider only the case < o < 2. 

Here, we will provide a systematic and controlled approach to the quantitative compu- 
tation of the thermodynamic momenta, usually used in numerical analysis. These momenta 
are related to the Binder's cumulant and to various thermodynamic functions like the sus- 
ceptibility. We will concentrate on the scaling properties of the coupling constants defining 
the system in the vicinity of the critical point. Our method is quite general and should 
apply to a large extend to the investigation of finite-size scaling in systems with long-range 
interaction in the vicinity of the critical point. 

The plan of the paper is as follows. In Section || we review, briefly, the (/9 4 -model with 
long-range interaction and discuss its bulk critical behavior. Section [H| is devoted to the 
explanation of the methods used here to achieve our analysis. We end the section with the 
computation of some thermodynamic quantities of interest. In Section IV we discuss our 
results briefly. In the remainder of the paper we present some details of the calculations of 
some formula used throughout the paper. 



II. FINITE-SIZE SCALING FOR SYSTEMS WITH LONG-RANGE 

INTERACTIONS 



In the vicinity of its critical point the Heisenberg model, with long-range interaction 



decaying as power-law, is equivalent to the d- dimensional (9(ra)-symmetric model |H| 



d d x 



(WV) 



2 2 1 4 



(2.1) 



where <p is a short hand notation for the space dependent n-component field f(x), r = r 0c +to 
(to oc T — T c ) and uo are model constants. V is the volume of the system. In equation (j2.1|), 
we assumed % = ks = 1 and the size scale is measured in units in which the velocity 
of excitations c = 1. We note that the first term in the model denotes k a \if(k)\ 2 in the 
momentum representation where the parameter < a < 2 takes into account short-range as 
well as long-range interactions. (3 is the inverse temperature. The nature of the spectrum 
suggests that the critical exponent r] = 2 — a ||IT| , |T2"|| . Here we will consider periodic boundary 
conditions. This means 
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<P(x) = -^7 eX P ( ife ' X ) ' ( 2 - 2 ) 



where fc is a discrete vector with components fcj = 27TOj/L (n$ = 0, ±1, ±2, • ■ ■ , i = 1, ■ • • , d) 
and a cutoff A ~ a -1 (a is the lattice spacing). In this paper, we are interested in the 
continuum limit i.e. a — > 0. As long as the system is finite we have to take into account the 
following assumptions L/a — ► oo, £5 — > 00 while is finite. 



Fisher et al. [1 1] and Yamazaki et al. [12|] have shown that for the model under consider- 
ation the Landau theory holds for d > 2cr. In the opposite case i.e. d < 2a an expansion in 
powers of e = 2a — d — 4 — d — 2r] takes place, where 2a plays the role of the upper critical 
dimension. We will present the renormalized parameters which characterize the bulk critical 
behavior and appear in the scaling functions. Since the computations are standard Q, we 
will be quite brief. 

The application of the renormalized theory, above the critical temperature, to the model 
Hamiltonian requires a scaling field amplitude Z, a coupling constant renormalization Z g 
and a renormalization of the p 2 insertions in the critical theory Z t . In term of these, we 
define as usual 

t = ZZ~\r -r 0c ) and g = £' £ Z 2 Z g 1 u . (2.3) 

In the remainder we will work in units where the reference length i is set to unity. To one 
loop order the renormalization constants in the minimal subtraction scheme are given by [12j 



Z = 1 + 0{g 2 ) (2.4a) 



Z t = l + —-g + 0(g 2 ) (2.4b) 



Z g = 1 + —g + 0(g 2 ) (2.4c) 



In Eqs.(pj) 



where ip(x) is the digamma function. 

The fixed point of the (3 function is at g = g* with 

Before starting to investigate the finite-size scaling in the field theoretical model under 
consideration, we shall recall briefly the corresponding renormalization group formalism. 
In the continuum limit, the lattice spacing completely disappears. The integration over 
wave vectors of the fluctuations are evaluated without cutoff and are convergent. When 



4 



some dimensions of the system are finite the integrals over the corresponding momenta are 
transformed into sums. Since the lattice spacing is taken to be zero, the limits of the sums 
still tend to infinity. 

From general renormalization group considerations an observable X, the susceptibility 
for example, will scale like ||15||: 



X[t, g, £, L] = ((p)X [t(p),g(p), £p, L] , (2.7) 

where t is the reduced temperature, g a dimensionless coupling constant and L the finite-size 
scale. The length scale £ is introduced in order to control the renormalization procedure. 

It is known that in the bulk limit, when g(p) approaches its stable fixed point g* then 
we have 

t(p) « tp 1 ^ and C(P) ~ P lx/ \ (2-8) 

where 7^ and v are the bulk critical exponents measuring the divergence of the observable 
X and the correlation length, respectively, in the vicinity of the critical point and p is a 
scaling parameter. Using dimensional analysis together with equation ( [2.7D one gets 



X[t,g,£,L] = ({p)X [t{p){p£)\g{p)^L/£p\ . (2.9) 

Choosing the arbitrary parameter p = L/£, we obtain the well known finite-size scaling 
result 

X[t,g,£,L] = L^f(tL 1 ^). (2.10) 

Here the function f(x) is a universal function of its argument. In the remainder of this 
paper we will verify the scaling relation ( p. 10 ) in the framework of model (|2.1|) . 



III. FINITE-SIZE SCALING BELOW THE UPPER CRITICAL DIMENSION 

A. Method 

The method, we shall use here to analyze the finite-size scaling of the model under 
consideration, is originally due to Luscher [|TBJ in his study on the quantum 0{n) nonlinear 
a model in 1 + 1 dimensions. An extension of the method was employed by Brezin and 
Zinn- Justin |T7J and by Rudnick, Guo, and Jasnow \TE\ in their works on the finite-size 



scaling in systems with short-range potentials. Very recently it was used in the investigation 
of crossovers in quantum 0(n) systems near their upper critical dimension [19|]. We will 



see here that the problem related to finite-size scaling in systems with long-range forces 
can be successfully analyzed by the same approach. Nevertheless, here we will observe the 
emergence of some subtleties, which need to be discussed. 

The central idea of the method is that at finite linear size L of the system, one can treat 
the k = mode of the field <p(x), playing the role of the magnetization, separately from the 
non zero k modes. The non-zero modes are treated perturbatively using the loop expansion. 
They are integrated out to yield an effective Hamiltonian for the lowest mode only. All 
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the modes being integrated out are regulated in the infrared by |fe| CT and consequently the 
process is necessarily free of infrared divergences. On the other hand the renormalizations of 
the bulk theory control the ultraviolet divergences at finite size. In other words if we define 
by 



0=i / d d xip(x) (3.1) 
V Jv 

the total spin by unit volume, then, from 7i, we can get an effective Hamiltonian function 
of 4> after entirely integrating out the (p(k ^ 0) fields: 

H cS = y ^ + ^0 4 ) . (3.2) 

The coupling constants R and U are computed in powers of e, with the initial coupling 
constants renormalized as in their bulk critical theory. This approach will rule out all the 
ultraviolet divergences of the bulk critical point. The new coupling constants are necessarily 
free of all ultraviolet divergences since the theory is superrenormalizable [12| . They are also 



free of infrared divergences as we are only integrating out finite modes. Obviously, these 
constants must obey the scaling forms, 

R = L^ 2 f R (tL l ' v ) and U = L d ~ 4+2f ? f v (tL 1 ^) (3.3) 



for t ^ 0, where fu and fu are scaling functions which are properties of the bulk critical 
point. They are analytic at t — 0. This is a consequence of the fact that only finite modes 
have been integrated out. 

Once the scaling functions fu and fu are known one can attack the problem of computing 
observables in the (p 4 theory with the action 7i e s- This theory is in dimension d close to the 
upper critical dimension 2a (not in d close to the usual 4), and the problem seems to be 
unsolvable. In the next section we will show that it is not the case. 

In order to investigate the long distance physics of the finite system, one has to calculate 
thermal averages with respect to the new effective Hamiltonian defined in ( |3.2|) . They are 
related to the thermodynamic functions of the system under consideration. The averages of 
the field 4> are defined by 

*»=(fr , )v J fif7 ( c T< r ) - (3.4) 

Using an appropriate rescaling of the field <f>: $ = ( UL J 0, we can transform the effective 
Hamiltonian into 

H cS = \z<$> 2 + U\ (3.5) 

where the scaling variable z = RL d l 2 XJ~ x l 2 is an important quantity in the investigations of 
finite-size scaling in critical statics [fT~T|JT~8[1 as well as in critical dynamics |2l1| , |2T|j . With the 
effective Hamiltonian (|3.5|), we obtain the general scaling relation 
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T P (d-4+2 V )/2 ( j{d-A+2r,)/2\ 

M2p = L-^+"> L m f 2r {Rt'-'.L——} (3.6) 

for the momenta of the field (p. Having in mind Eqs. (|3.3| ), we can write down Eq. (|3.6|) in 
the following scaling form 

M 2p = L-^-^F^tL 1 /"), (3.7) 

in agreement with the finite-size scaling predictions of ( 2.10| ). In eq. ( |3.7| ), the function 
T2 P {x) are universal. 

All the measurable thermodynamic quantities can be obtained from the momenta M.2p- 
For example the susceptibility is obtained from 



X 



-[ d d x (<p(x)(p(0)) = L 2 -vT 2 {tl}l v ). (3.8) 
n Jv 



Another quantity of importance for numerical analysis of the finite-size scaling theory is the 
Binder's cumulant defined by 



1M 



B = l -3M (3 ' 9) 

In the remainder of this section we concentrate on the computation of the coupling con- 
stants R and U of the effective Hamiltonian (|3.2|) for the system with long-range interaction 
decaying with the distance as a power law. As a consequence we will deduce results for 
the characteristic variable z = RU^ 1 ! 2 ' L 2 ~' t] ~ e l 2 , the susceptibility x an d the amplitude ratio 
r = Ai^/Ai 2 . entering the definition of the Binder's cumulant. 



B. Computation of the coupling constants R and U 

As we explained above, loop corrections will be treated perturbatively on the non-zero k 
modes. At the tree level (lowest order in e = 2a — d) this procedure generates a shift of the 
critical temperature T c and a change of the coupling constant u and additional operators 
involving powers of if larger than 4. The calculations will be performed in the renormalized 
theory. The renormalized coupling constant ur is expressed in terms of the dimensionless 
coupling constant g = 1 £ ur in which the parameter i is an arbitrary length scale. Here 
we will work in system in which I = 1. Throughout these calculations we use the minimal 
subtraction scheme. In this scheme, the counterterms of the massless theory including the 
if 2 insertions are introduced. The one-loop counterterm for the coupling constant and the 
ip 2 insertion will be the only one relevant in the lowest corrections. 

The finite-size correction to the renormalized coupling constant t is given by 

K.{t,g,L) = (^ + 2)^E / 7^7 ( 3 - 10 ) 

to one-loop order. 
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In order to investigate the finite-size scaling of the model under consideration one can 
use a suitable approach allowing to simplify the analytical calculations. In the case o = 2 it 
is possible to replace the summand by its Laplace transform. This is the so called Schwinger 
representation. The aim of this approach is to reduce the d- dimensional sum in the r.h.s of 
equation ( |3.10| ) to the one-dimensional effective problem. In the general case of arbitrary 
a, one cannot just use the Schwinger transformation or at least in its familiar form. So we 
have to solve the problem by introducing some kind of generalization for it. In the spirit 
of the same problem a method to investigate the finite-size scaling in the framework of the 
mean spherical model was suggested in Ref. [021. The method is based upon the following 
genius identity 



= / dxexp (— xz) x a ~ 1 E aa (— x a ) 

1 + z a Jo 



where the functions 



(3.11a) 



(3.11b) 



is the so called Mittag-Leffler type functions. For a more recent review on these functions 
and other related to them, and their application in statistical and continuum mechanics see 
Ref. [|23J. See also Ref. and Appendix 0. 

Using the identity (]3.11|) , one gets, after some algebra, 



WL(t, g, L) — (n + 2)g—^-_ / dxx^E*.* 



(2vr) CT Jo 



^/2 



(2tt)' 



A\x) 



(3.12a) 



where 



A{x) 



£=-oo 



d 2 



(3.12b) 



The analytic properties of the function A(x) are known very well. For large x, A(x) — 1 
decreases exponentially and the integral in the r.h.s of Eq. ( p. 12a ) converges at infinity. For 

small x, the Poisson transformation A(x) = (^j 2 A shows that A(x) converges. 

For small x the integral in the r.h.s of Eq. (|3.12a|) has ultra violet divergence for Re d > a. 
So, an analytic continuation in d is required to give a meaning to the integral. Adding and 
subtracting the small asymptotic behavior of the function A(x), we get after some algebra 



" cr—d 



WjJt, g, L) = (n + 2) g7 —-F d , a (tL° 



(27r) ff 
+2Tr(n + 2)gL a - d 



i -i 



{tL c 



id/o--l 



where 



F d ,a (y) 



dxx 2 



--1 



2 ' 2 



yx 



a/2 ' 



(2tt) ct 



A d (x) - 1 



x) 



(3.13a) 



(3.13b) 
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In the particular case a = 2, from Eq. (|3.13|) we recover the result of Ref. |T7 . 

By introducing the (f 2 counterterm insertion the renormalized coupling constant t is 
replaced by tZ t , where Z t is given by ( |2.4b| ). Hence to one loop order we have 

n + 2" 



R = t(l+g 



+ W t d Jt,g,L). 



(3.14) 



At d = 2a, W d a (t, g, L) has a simple pole. An expansion about this pole leads to the 
final expression 



R = t + 



n + 2 



a 



gtlnt + 2 a ~\n + 2)T(a)gL-° "F 2(TiCT (tL a ) + O [g 



(3.15) 



This result shows that, at the critical point, R has the required scaling properties of Eq. Q3.3|), 
since 



i n + 2 „. 
y- 1 = a —e + 0{e 2 ). 



n + 8 



(3.16) 



For the finite system the renormalized coupling constant g, to one-loop order, is shifted 
by a quantity expressed in the form 



(3.17) 



As one can see the summand here can be expressed as the first derivative of the summand 
of Eq. (|3.10|) with respect to t. So, the result for U can be derived from that of R. Using 
this fact one gets 



W 9 d Jt,g,L) = (n + 8)g 2 



(27T) 



a(47r) d / 2 



r(d/2) 



(3.18) 



where the prime indicates that we have the derivative of the function F with respect to its 
argument. 

At the fixed point one ends up with 



U = g 



I 9— (1 + ^ t) + 9^T(a)Fi a>a (tL°) + O (g 2 ) 



a 



(3.19) 



for the renormalized coupling constant U. Eq. (|3.19|) is obtained using the fact that at 
one-loop order the coupling constant is renormalized by Z g form Eq. ( |2.4c ). The obtained 
expression ( p,19|) shows that the coupling constant U obeys the scaling law of Eq. (p.3|). 
Note that U has a finite limit as t — > 0, i.e. it is analytic at the bulk critical temperature. 
Indeed as t — > one can use the expansion of the function F2 a ,a(y) for small y given by (see 
Appendix [5]) 



F 2 aAv) =F 2 ^(0) + 2- a yC a - 



)i-<j 



aT(a) 



ylny + 0(y 2 



(3.20a) 
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where 



C„ 



du 



Y(a) Jo u 



Es. i 

9 ! 



u 



a/2 



(2tiY 



—A 2 °{u) + — 



(3.20b) 



After substitution of ( p. 20] ) in ( |3.19|) the terms proportional to logy cancel, which shows 
that the coupling constant U is finite at t = 0. Whence, one gets 



U = gL- 



1+^^ (i + fl»C (T ) Our 



showing that £7 is analytic, as it should be, at the critical point. 



C. Some thermodynamic quantities 



1. Shift of the critical point 



It is obvious that the coupling constant R in the effective Hamiltonian ( |3.2|) is just the 
deviation of the temperature of the system from its 'critical' value. By setting t = in 
( p,15| ), we obtain an expression for the finite-size shift of the bulk critical temperature T c . 
This is given by 



n + 2 

T c - T C {L) = £ 2°- 1 ——r(a)L-°F 2 ^(0), 
n + 8 



(3.21) 



where the coefficient F 2<T)CT (0), appearing in the right hand side of ( |3.21| ) can be evaluated 



for some particular values of the interparticle interaction range a [24 

r 2c (1/2) , 

4C (1/2)/? (1/2), 
-4.82271993, 
-8 In 2, 



F 2 *,a (0) 



a = 1/2, 
a = l, 
a = 3/2, 
a = 2. 



(3.22) 



Here ((x) is the Riemann zeta function with ( = —1.460354508... and (3(x) is the analytic 
continuation of the Dirichlet series: 

m= ]^ (2e + iy> 

with /3 (|) = 0.667691457.... Remark that the function F 2(Tj(T (0) increases as the parameter 
a vanishes. 

In fact, since there is no true phase transition in the finite system, the critical temperature 
is shifted to a 'pseudocritical' temperature, T C (L), corresponding to the rounding of the 
thermodynamic singularities holding in the bulk limit. From fl3.21| ) one remarks that T C (L) 
is larger than T c , confirming previously obtained results in the framework of the spherical 



model p5|. Notice also that for the shift exponent A, we get A = a to lowest order in e. 
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2. Binder's cumulant 



In this subsection we are interested in the calculation of the amplitude ratio r = Ai^/Ai 2 , 
instead of the Binder's cumulant from definition (|3.9|) . This quantity can be expressed in 
power series of the scaling variable z = RL 2 ~ r, ~ £ l 2 U~ l l 2 as 



n 



(H 



4F» (i(n + 2] 

r(i(n + 6))r(| 



+z 



r(|(n + 6)) T(\(n + 2)) r(j(n + 4) 



r(|(n + 4)) r(|n) 
(J(" + 4 ) 



2) 



+ 3- 



r(j(n + 4))r(jn) n(l(n + 2 



n 




(3.23) 



So, in order to obtain a result for r it is enough to evaluate z at the fixed point g* and to 
deduce the value for the Binder's cumulant. As we mentioned before this parameter appear 
in all thermodynamic functions through the momenta defined earlier in this paper. 
At the fixed point g* in the vicinity of the upper critical dimension, we obtain 



RL 2 -v 



VUL 
1 



'9 



fixedpoint 

la 



n - 4 



n + 8 



In y) +2 a - 1 e 



n + 2 
n + 8 



-e2°- 2 yT(a)F> (y) 



(3.24) 



This result is obtained by using ( |3.16[) and the fact that up to one loop order the terms 
proportional to InL cancel. In Eq. ( |3.24| ), we introduce the scaling variable y = tL 1 ^. 
Finally let us notice that from this equation one can see easily that z* verifies the finite-size 
scaling hypotheses and consequently all the thermodynamic functions do. 
At the critical temperature T c (i.e. t — 0, and so y = 0), we obtain 



n + 2 




iW(0) + 0{e) 



(3.25) 



Numerical values for the amplitude ratio (|3.23| ) can be obtained by replacing the value 
of Zq form ( [3.25 ) and taking some specific values of the small parameter e. Note that the 
scaling variable z is proportional to y/e as it was found previously (see Ref. [IT] for example) 
in the case of short-range forces. Furthermore it coincides with the result of Ref. for the 
scaling variable x in the case of long-range interaction. Consequently all the thermodynamic 
function will be computed in powers of y/e. 



3. Magnetic Susceptibility 

As we mentioned earlier, there is no phase transition in the finite system under consid- 
eration. Consequently there will be no 'true' correlation length. An expression for it can be 
deduced from that of the susceptibility fl3.8Q trough the relation: 
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2-7/ 



(3.26) 



The analyticity of the susceptibility is a consequence of that the coupling constants R and 
17. 

From ( |3.8|) in the region tL a <C 1 (i.e. 2 1), we obtain for the susceptibility 



X 



2V2 ^+%T(\{n + 2 
r (a) 



y / (47r)-r(a) ra 

l_ n T 2 (i(n + 2 



l-z — 



r(f) 



.n + 8 



4 



r 2 (f) 



4 T (\(n + 2\ 



r(f) 



i + ^i>)c) +o(^,^ 



(3.27) 



at the bulk critical point T c . 

To the lowest order in e, after taking the limit n — > 00, we find that the correlation 
length scales like 

t _-l/2<rr 



confirming the results obtained in the spherical model |26| and showing that this behavior 
is not a characteristic of the spherical limit i.e. n — > 00. 

In the region tU 1 (i.e. z 3> 1), from Eq. ( |3.8| ), we get 



X 



1 

t 

1 

— [n 
2 



n + 2 



(7 



£ In t - 2°-\n + 2)r(<r)$ {tL^ 1 F 2a>a (tU) 



The function F d>(T (y) has the following large y asymptotic behavior (see Appendix 



^3.28) 



F d Ay) - 

for the case < o < 2, and 

^,2(2/) 



?/ 2 r(-|) 



|Z| d+<J 



(3.29a) 



4vr 2 



+ d(27r)( 5 - d )/ 2 2 /^ 3 )/ 4 e-v / ^ 



(3.29b) 



for the particular case a = 2. These results show that the last term in Eq. ( p.28|) is just 
canceled by the first term in Eqs. ( |3.29| ). 

In the case of long-range interaction < a < 2, we obtain for the susceptibility 



X = Xc 



1 - ag(n + 2)2 3 - 2 {tL^ ^^ E' ^ + Off 



(3.30) 



in agreement with the finite-size scaling hypothesis ( |3.8|) . Eq. (|3.30|) shows that the finite- 
size scaling behavior of the system is dominated by the bulk critical behavior, with small 
correction in powers of L. It should be noted that the above result cannot be continued 
smoothly to the case of short-range interaction o = 2, since then F^ 2 {y) (see Eq. ( |3.29b| )) 
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falls off exponentially fast and, correspondingly, the the finite-size corrections to x are ex_ 
ponentially small: 



X = Xc 



1 - 8gV27T(n 



+ 2) (tL 2 ) 



(3.31) 



At this point we are in disagreement with the statement given in Ref. [27| that the approach 
used in Refs. |PJT7|-pT| yields an incorrect non-exponential result. Note that all our calcula- 
tions are up to the order g l . It is interesting to see what happen in higher order, e.g. g 2 , in 
this case, however, we need to have at our disposal the corresponding high order terms in 
Eqs. ( |3.15D and ( |3.19|) . Indeed it is beyond the scope of the present study. First the power 
law fall off of the finite-size corrections to the bulk critical behavior, due to long-range na- 
ture of the interaction, was found in the framework of the spherical model p8| , p9| . Here 
extended this result to finite n using a perturbative approach. 



we 



IV. CONCLUSIONS 

In this paper, we have investigated the finite-size scaling properties in the 0(n)- 
symmetric if 4 model with long-range interaction potential decaying algebraically with the 
interparticle distance. We have found that the methods developed in Refs. fl6Hl9f can be 
successfully extended to systems with long-range interaction by combining them with other 
known techniques. These techniques allow the investigation to be simplified and express the 
results for various thermodynamic functions in terms of simple and known mathematical 
functions. 

Here we restricted our calculations to the critical domain T > T c and investigated the 
model in dimensions less than the upper critical one, which turns out to be 2a (0 < a < 2). 
We constructed an effective Hamiltonian, from the initial one, with new coupling constants R 
and U. These constants obey the scaling hypothesis (|3.3|). We found that the even momenta 
of the field ip, related to the thermodynamics of the finite system, are scaling functions of 
the characteristic variable 

z = RU- l ' 2 L 2 -'^' 2 . 

This variable has the required scaling form predicted by the finite-size scaling theory. From 
the obtained forms of the constant R and U one concludes that z is a universal quantity, 
which does not depend of the details of the model. 

We evaluated the finite-size shift, the susceptibility and the amplitude ratio r = M. A / ' }A\ 
at the tree level (lowest order in e). We observed that the critical behavior of the system is 
dominated by its bulk critical behavior away from the critical domain and that the finite-size 
scaling is relevant in the vicinity of the critical point. The amplitude ratio r is evaluated 
as an expansion in powers of z ~ y/e. Our result is in consistency with that of Ref. ||. 
But it is disagreement with the numerical results of the same paper. There, it has been 
found, using the Monte Carlo method, that r has an expansion in e instead of its square 
root. At this time, we do not have a reasonable explanation of this fact. It is also possible 
that higher order in e could improve the result. An amelioration of the result could also 
come from accounting finite cutoff effects, which were to be relevant in the investigation of 
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finite systems and the comparison of the results with numerical works |30],|3T| . However this 
is the subject of another publication. 

Notice that in the only work devoted to the exploration of finite-size scaling in 0{n) sys- 
tems with long-range interaction (Ref. |J) the pertinent integrals have to be evaluated only 
numerically, due to the choice of a parametrization that does not reduce the <i-dimesional 
problem to the effective one dimensional one. The approach we used here is more efficient 
in the sense that the corresponding final expressions can be handled by analytical means. 
Consequently, we cannot make a direct comparison between the results of this paper and 
those obtained there. 

Let us note that it would be interesting and useful to extend the result obtained here in 
the static limit to models including dynamics, since we believe that this is closely related to 
the extensively investigated filed of quantum critical points i.e. phase transitions occurring 
at zero-temperature. In particular we find it useful to investigate the critical dynamic of the 
quantum model considered in Ref JSJ in the large n limit. 
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APPENDIX A: SOME PROPERTIES OF THE MITTAG-LEFFLER TYPE 

FUNCTIONS 

The Mittag-leffler type functions are defined by the power series p3| : 

oo z k 

KAz) ^f^kWY a ' p>0 - (A1) 

They are entire functions of finite order of growth. Let us mention that the function cor- 
responding the particular case j3 = 1 was introduced by Mittag-Leffler. These function are 
very popular in the field of fractional calculus (for a recent review see Ref. [fffil ). 

One of the most striking properties of these functions is that they obey the following 
useful identity |23| ) 



1 



dxe~ x x p - l E a ^ (-x a z) , (A2) 



l + z 

which is obtained by means of term- by-term integration of the series ( |A1|) . The integral in 



Eq. (|A2|) converges in the complex plane to the left of the line Kez = 1 1//q , | argz| < ^an. 
The identity ( |A^ ) lies in the basis of the mathematical investigation of finite-size scaling 
in the spherical model with algebraically decaying long-range interaction (see Ref. @ and 
references therein). 

In some particular cases the functions E a ^(z) reduces to known functions. For example, 
in the case corresponding to the short range case we have: 
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#i,i(z) = exp(z). (A3) 
Setting z = y~ a , y > 0, and x = ty, we obtain the Laplace transform 



1 + Z Q 



dte- zt t p - l E a) p {-t a 



(A4) 



from which we derive the identity ( |3.11| ) by setting (3 = a. 

The asymptotic behavior of the Mittag-LefHer functions is given by the Lemma p3[ 



Let < a < 2, j3 be an arbitrary complex number, and 7 be a real number obeying the 
condition 

1 r -I 

-air < 7 < mm{7r, air}. 
2 

Then for any integer p > 1 the following asymptotic expressions hold when |z| — ► 00: 
• At I argz| < 7, 



a 



1 ~(1-/3)/"^ 1/q " ; 



£ r(/3 - a*) 



-p-i 



(A5) 



At 7 < I argz| < 71, 



J r(/3 - afc) 



O \z 



-p-i 



(A6) 



APPENDIX B: ASYMPTOTIC BEHAVIOR OF THE FUNCTION F d)G {y) 
To obtain the small y behavior (|3.20a| ) of the function Fd, a (y) we use the identity [^2 



roo Arr 

\ R< f > = a — [E aA (-x a ) - E aA (-yx a )] 
Jo X 



and the definition of the function F dj(T (y): 



poo 

F d ,a (y) = / dxx^E 
Jo 



yx 



a/2 ' 



2 > 2 



(2tt) <j / 



^fx) - 1 



n y/2 
x) 



After some algebra one obtains: 



F 2 aAv) =F 2 ^{Q) + 2-°yC a 



)1-<J 



-yln^ + Od/ 2 



where 



T{a) Jo u 



Es. 1 



ct/2 



(27r) CT / TT' 7 



— A 2a (u) + — 



(Bl) 



(B2) 



(B3a) 



(B3b) 



15 



To obtain the large y asymptotic behavior ( |3.29| ) of the function Fd >(T (y) we rewrite ( p2| 
in the form 



' y ' Jo 2,2 v (27r) <T yY ^ 2,2 V ( 2? 



cr/2 v 

(2^ 



(B4) 

Using the identity 

poo a / / \ 

y dxx^E'j (-x a/2 ) =1, cr > (B5) 



From the second term of Eq. flB4|) we obtain the first terms of Eqs. ( p.29a ) and (|3.29b[) 
respectively. 

Next taking into account Eq. (|A6|) or Eq. (|A3|) for the function E a ^(z) and after 
subsequent integration in the first term of Eq. (|B4j), we obtain finally the asymptotic 
behavior given by Eqs. ( |3.29aj) and (|3.29b| ). 
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